FERMIL AB-PUB- 1 2-466-T , ANL-HEP-PR-12-58, EFI-12-20, MAN/HEP/2012/11 
KCL-PH-TH/2012-30, LCTS/2012-15, CERN-PH-TH/2012-195, CNU-HEP-12-01 



(N 

o 

(N 



(N 
(N 



CPsuperH2.3: 

an Updated Tool for Phenomenology in the MSSM 
with Explicit CP Violation 



q\ J. S. Lee a ' 6 , M. Carena c , J. Ellis d ' e , A. Pilaftsis / and C. E. M. Wagner 9 '' 1 

CN . a Department of Physics, Chonnam National University, 

300 Yongbong-dong, Buk-gu, Gwangju, 500-151, Republic of Korea 

43 ■ 

o ; Department of Physics, National Tsing Hua University, Hsinchu, Taiwan 300 

c Fermilab, P.O. Box 500, Batavia IL 60510, U.S.A. 

d Theoretical Particle Physics and Cosmology Group, Department of Physics, 
King's College London, London WC2R 2LS, United Kingdom 

(N ; 

£> 1 e Theory Division, CERN, CH-1211 Geneva 23, Switzerland 



f Consortium for Fundamental Physics, School of Physics and Astronomy 
University of Manchester, Manchester M13 9PL, United Kingdom 



9 HEP Division, Argonne National Laboratory, 9100 Cass Ave., Argonne, IL 60439, USA 

O 

CN . Enrico Fermi Institute, Univ. of Chicago, 5640 Ellis Ave., Chicago, IL 60631, USA 

^ ; ABSTRACT 

We describe the Fortran code CPsuperH2 . 3, which incorporates the following updates compared 
with its predecessor CPsuperH2 . 0. It implements improved calculations of the Higgs-boson masses 
and mixing including stau contributions and finite threshold effects on the tau-lepton Yukawa 
coupling. It incorporates the LEP limits on the processes e + e~ — > HiZ, HiHj and the CMS limits 
on Hi — > tt obtained from 4.6 fb _1 of data at a centre-of-mass energy of 7 TeV. It also includes 
the decay mode Hi — > Zj and the Schiff-moment contributions to the electric dipole moments 
of Mercury and Radium 225, with several calculational options for the case of Mercury. These 
additions make CPsuperH2 . 3 a suitable tool for analyzing possible CP-violating effects in the 
MSSM in the era of the LHC and a new generation of EDM experiments Q. 
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1 Introduction 



Supersymmetry is one of the most attractive possible scenarios for physics beyond the 
Standard Model in the TeV energy range, and the minimal supersymmetric extension of 
the Standard Model (MSSM) is the simplest framework that incorporates this scenario. 
The MSSM allows for many possible CP- violating phases, notably those in the SU(3), 
SU(2) and U(l) gaugino masses M^i and in the trilinear soft supersymmetry-breaking 
parameters associated with the third-generation Yukawa couplings At 7 b, T - 

CPsuperH jTJ[2] is an evolving software tool that incorporates the effects of these CP- 
violating phases into the calculation of Higgs boson masses, couplings and other properties. 
The previous version, CPsuperH2 . |T] , also incorporated a number of 5-physics observ- 
ables, including the branching ratios for B s — > Bd — > t + t~, B u — > ru, B — > X s ~f 

and the CP- violating asymmetry in the latter decay, Acp, as well as the supersymmetric 
contributions to the B® d — B® d mass differences. Also included in CPsuperH2.0 were cal- 
culations of two-loop supersymmetric Higgs contributions to the electric dipole moments 
(EDMs) of Thallium, the electron and muon Q. 

After the recent discovery at the LHC [1] of a new particle resembling a Standard- 
Model-like Higgs boson, Higgs physics and searches for associated CP-violating effects are 
entering a new era, in which the LHC experiments are now also probing directly the possible 
existences of heavier supersymmetric Higgs bosons. In parallel, the LHCb experiment is 
taking experiments on 5-physics to a new level of precision, and a new round of experiments 
is set to improve significantly sensitivities to a wider range of EDMs, including also Mercury 
and Radium. This juncture in the exploration of physics at the TeV scale, both direct 
and indirect, is the appropriate moment to document the capabilities of a new update of 
CPsuperH. 

The followings are the main new features incorporated in CPsuperH2 . 3. As described 
in Section 2, it features an improved treatment of CP-violating effects on Higgs-boson 
masses and mixing, with stau contributions included as in Ref. [5]. We also take consistently 
into account finite radiative effects on the tau-lepton Yukawa coupling. As described in 
Section 3, CPsuperH2 . 3 also features an implementation of the LEP limits on the processes 
e + e~ — > HiZ,HiHj [6], as well as the limits on Hi — > fr obtained by CMS with 4.6 fb _1 
of LHC data at a centre-of-mass energy of 7 TeV [7J, where Hi (with i = 1,2,3) denotes 
collectively the three neutral Higgs bosons #1,2,3 in the CP-violating MSSM. In the same 
section, we also present theoretical predictions for Hi — > Zj, as these decay modes can 
become detectable as the integrated LHC luminosity increases. Section 4 outlines how 
CPsuperH2.3 incorporates calculations of the EDMs of Mercury and 225 Ra that include 
estimates of the contributions due to Schiff moments [5]. Each section includes figures that 
illustrate some typical results obtained using CPsuperH2.3. As explained in Section 5, 
an SLHA2 interface was created to facilitate the comparison and linkage of output data 
between CPsuperH2 . 3 and other public codes, in accordance with the SUSY Les Houches 
Accord [T0j[TT]. Finally, Section 6 summarizes our updates to the code. 

'For an alternative code treating similar physics, see [3]. 
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2 Stau Contributions to Higgs masses and Mixing 



The scalar tau contributions to the masses and mixing of the neutral Higgs bosons have 
been included as in Ref. [5], similarly to the sbottom corrections with 
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More precisely, in the (0i,0 2 ,o) T basis, the scalar tau contributions to the entries of the 
neutral Higgs mass squared matrix from the one-loop effective potential are given by 
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The tan /3-enhanced terms containing /i and v4 r are included in R T and given by 

|/x| 2 tan/3 — 3?e(A T /i) , |y4 r | 2 — 3f?e(A T/ u) tan/3 



and the loop function 
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■^To match the CPsuperH convention, the sign of A T is flipped compared to Ref. [5]- 
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Since the stau contributions to the mass splitting M\± — M\ are small, we have neglected 
them. For the same reason, we set 
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Note that in the limit of small mass splitting one has 
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with M| = [m\ + m 2 v )/2. So, in this limit, the leading contributions proportional to tan 4 /3 
are given by 
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which are the same as those presented in [T^HT^] . Note that here we take into account the 
fi j2 mass splitting. 

In Fig. [H we demonstrate the impact of the corrections from the light scalar taus for 
large /i = 1 TeV for several values of tan/3. When m~ ~ 1 TeV, the stau corrections 
are negligible and the lightest Higgs boson mass increases by the small amount of ~ 0.05 
GeV even when tan/3 = 60. However, these corrections are larger for lighter staus, and 
the lightest Higgs boson mass may decrease by as much as ~ 4 GeV, if tan/3 = 60 and 
m~ ~ 100 GeV. For comparison, as shown in Table IU the masses and mixing matrix of 
the neutral Higgs bosons without including the stau effects are stored in the array RAUX_H: 



RAUX_H(511 — 513) : Three Higgs masses without including the stau effects; 

RAUX_H(520 — 528) : The nine elements of the mixing matrix O a i without including 
the stau effects. 



3 Collider Limits 



The limits on Higgs boson production at LEP and the LHC experiments have been imple- 
mented as follows. 



4 



128.5 



MHMAX scenario but with (j, = 1 TeV and varying M L3 =M E3 



tan(3 = 40 



123.5 




M„, = 1 TeV 



100 200 300 400 500 600 700 800 900 1000 



m_ Tl [GeV] 



Figure 1: The lightest Higgs boson mass as a function of the lighter stau mass for several 
values o/tan/3: tan/3 = 10 (black), tan/3 = 20 (red), tan/3 = 30 (green), tan/3 = 40 (blue), 
tan (3 = 50 (yellow), and tan /3 = 60 (magenta) , as obtained varying = m_g 3 wift /ixed 
fi = 1 TeV. We use the MHMAX scenario for the other supersymmetric model parameters: 
m Qi = m(Ji = mt)i = Msusy = 1 TeV, M 1 = 100 GeV , M 2 = 200 GeV , M 3 = 800 GeV, 
and A t = \/6 M SU sy + /V tan/3 with A b = A T = A t . We have fixed M H ± = 1 TeV. The thin 
short line around m~ = 1 TeV for each tan /3 value shows the lightest Higgs boson mass 
before the inclusion of the stau effects. 

3.1 LEP Limits 

We have implemented the LEP limits from the processes e + e~ — > HiZ and e + e~ — > HiHj 
using Table 14 (a) and Table 17 (a) of Ref. [6], respectively, assuming that the Higgs decay 
patterns are not drastically different from those in the Standard Model and in the MHMAX 
scenario with tan/3 = 10. We have required that each of the Higgs couplings to the gauge 
boson(s) normalized to the Standard Model values, gjj^v an d 9 < H i H :j zi should be smaller 
than the corresponding values in the Tables. The result is saved in 

- RAUX_H(430)=ILEP : (Excluded) or 1 (Allowed). 



For the more general cases, we refer to more refined tools such as HiggsBounds [15] 
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3.2 LHC Limits 



We have also incorporated the recent CMS limit from the search for the MSSM neutral 
Higgs bosons decaying into tau pairs based on 4.6 fb _1 [7j, and the result is saved in 

- RAUX_H(440)= ILHC7f 6 ^ rr : (Excluded) or 1 (Allowed) . 

In our implementation of the CMS limit, assuming the same K factors as in the SM, 
we approximate the production cross sections of the neutral Higgs bosons at the LHC as 
follows: 

/-pMSSMX LO 

a^( PP ^HiX) ~ I tZgsr a™(pp -+ H SU X) _ , (13) 



H SM H i 



where V = ggHi, bbHi, VVHi specifies each production process. The SM production cross 
sections are calculated by using HIGLU [16], BBHONNLO [IT] and HAWK [18] § In leading order, 
the process-dependent ratios are given by: 

— Gluon fusion 

ffl^ _ \S?(M Hi )f + \P?(M Hi )\* 

loJ 4 \s 9 SM (M Hl w - Hm " {) 

where S^ m (MhJ = J2f=b,t -^s/( T i/)- These factors are saved in the array CAUX_H (221-223) 
as shown in Table [6j 

— b- quark fusion 

2 

/pMSSM\ LO 2 Q P _ 

= R Hibb . (15) 



/-nMSSMX LO 




f 1 bbH t \ 




I F SM 1 


dHibb 


\ L bbHi J i 





where the extra factor in the second term accounts for the phase-space difference 
between scalar and pseudo-scalar decays. 



Vector-boson fusion 

/pMSSM\ LO 
I 1 VVHi \ 2 



. 1 VVHi } i 



9h,vv ( 16 ) 



In Figs. [2H31T41 we show the three process-dependent ratios as functions of the corresponding 
Higgs masses, assuming the MHMAX scenario. For completeness, also shown in Fig. [5] are the 
ratios of the Higgs-boson couplings to two photons defined by 

\s]{m h ^ + \p:\m h ^ 



^We thank Junya Nakamura for providing the SM cross sections. 
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Figure 2: The couplings of the MSSM Higgs bosons $ to two gluons normalized to the 
Standard Model values, PH igg P4\ ), functions of their masses for several values of tan (3. 
In the CP-conserving limit, $ = h, H, A with h(H) and A being the lighter (heavier) CP- 
even and CP-odd Higgs bosons, respectively. We assume the MHMAX scenario: m 



m 



u 3 



m 



d 3 



m 



L 3 ~ m E 3 

800 GeV, and A t 



200 GeV , M 3 



Q3 

Msusy = 1 TeV, fi = 200 GeV,Mi = 100 GeV , M 2 
V^Msusy + /V tan/3 with Af, = A T = A t . 



where S^ M (M Hi ) = 2J2f=b,t,c,r N c Qj F sf (T if ) - F sf (T iW ) } which are saved in the array 
CAUXJK231-233) as shown in Table EJ 

We observe that the branching ratio B(H\ — > 77) may be enhanced if the lighter 
stau has a low mass, particularly for large tan /3 as shown in Fig. |6j This enhancement is 
consistent with the signals observed [1], which may be larger than in the Standard Model. 
Furthermore, in Fig. [7] we show the correlation of B(H\ — > 77) with B(H\ — > Zj) in the 
MSSM, after the branching ratios are normalized to their SM predictions. Calculational 
details of B(Hi — > Zj) are described in Appendix A. 
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MHMAX scenario 
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Figure 3: The same as in Fig. [H but for the couplings of the Higgs bosons to b quarks fldV . 



Once the production cross sections of the neutral MSSM Higgs bosons have been 
calculated, we may require the sum 



a(pp -»• HiX) ■ B(Hi 



rr)+ £ 



a(pp ->■ <f>X) ■ B{ 



tt) 



Ma-M h .\<SM 



to be smaller than the observed limit for a given value Mn t of a specified neutral Higgs 
boson Hi. Here we are adding all Higgs production cross sections of <fi ^ Hi to the one of 
the Hi boson when their mass difference is smaller than SM. We take 5M = 0.21 * 130/2 ~ 
13 GeV, corresponding to the tau-pair mass resolution of ~ 21 % at a Higgs boson mass 
of 130 GeV [7]. Such a simple treatment is fairly accurate, as long as SM 3> T^, namely 
in the absence of strongly overlapping Higgs resonances [9]. Under these considerations, 
Fig. [8] shows the LEP- and LHC-excluded regions in the M4-tan /3 plane for the MHMAX 
scenario. We find our results reasonably consistent with those presented in Fig. 4 of [7J 
when M A > 90 GeV. 

Finally, for the value given by M HsM =SMPARA_H(20) , the decay widths and branching 
ratios of the SM Higgs have been calculated, and the results are stored in RAUX_H (600-700) , 



^In the CP-conserving limit, the CP-odd state A is Hi with 2 ai = 1. In the CP-violating case, we 
identify the "CP-odd state" as the one which has the the largest CP-odd component squared 2 ai . 
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MHMAX scenario 
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Figure 4: The same as in Fig. [H but for the couplings of the Higgs bosons to two vector 
bosons (Tb}) . 

see Table HI For comparison with other works, we display in Fig. [9] the decay widths into 
b quarks and photons (upper) and the branching ratios into W and Z bosons (lower) as 
functions of the SM Higgs-boson mass. We note that decay widths into off-shell vector 
bosons using the four-body phase-space have been implemented in CPsuperH2 . 3. 



4 Mercury and Radium Electric Dipole Moments 



CPsuperH2 . included two-loop Higgs-mediated contributions to the Thallium and electron 
EDMs [I], and CPsuperH2.3 extends these results to include calculations of the Mercury 
and 225 Ra EDMs that incorporate Schiff-moment contributions. In the case of the Mercury 
EDM, these are parameterized as follows 
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Figure 5: The same as in Fig. [H but for the couplings of the Higgs bosons to two photons. 

where d^g 1 ' 111 ' 1 [S] denote the following four different Schiff- moment induced Mercury EDM 
calculations: 

- LSxlO^e^/GeV, 
d£ g [S] ~ 7.6 x 10- 6 e^/GeV + 1.0x 10~ 3 e£ ] NN /GeV , 
d™[S] ~ 1.3 x 10- 4 e^/GeV + 1.4 x 10~ 3 e£l N /GeV , 
d™[S] ~ 3.1 x lO" 4 ^ ^ /GeV + 9.5 x 10~ 5 e£^ N /GeV, 

where the g^NN are the CP-odd irNN couplings, and we refer to Ref. [8] for further details. 
We note that is basically the same as that calculated in CPsuperH2.0. In the case of 
the Radium 225 Ra, we estimate the EDM through [8] 

d Ra ~ d Ra [S] ~ -8.7 x 10" 2 egf NN /GeV + 3.5 x 10" 1 /GeV . (19) 

Fig. [10] shows the four estimates of the Mercury EDM and the Radium EDM as func- 
tions of the parameter p that parameterizes the hierarchy between the first two and third 
generations in the CPX scenario. The new EDMs are stored in the array RAUX_H: 

- RAUX_H(411) = d^ g , RAUX_H(412) = d£ g , RAUX_H(413) = dg, RAUX_H(414) = 

10 
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Figure 6: The ratio of the MSSM-to-SM branching ratios of the lightest Higgs boson H\ 
decaying into photons as a function of the lighter stau mass. The scenario and lines are 
the same as in Fig. Ul 

- RAUX_H(420) = d Ra 
as shown in Table El 

5 SLHA2 interface 1 

CPsuperH2 . 3 also provides an output in accordance with the SUSY Les Houches Accords 
1 (SLHA1) [10] and 2 (SLHA2) [IT]. By taking 

• IFLAG_H(30) = 1 

in the run file, the output file cpsuper2 . 3_slha2 . out is generated. The output file of the 
current version includes the following blocks: 

"We thank Alexander Pukhov for helpful discussions regarding the interface. 
**Note that the default setting is IFLAG_H(30)=0, which is not recommended for generating files when 
a scan of the parameter space is performed. 
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Figure 7: Correlations of the MSSM-to-SM branching ratio of the lightest Higgs boson H± 
decaying into two photons with the respective one for the decay Hi into a Z boson and a 
photon, for discrete choices o/tan/3. The scenario and lines are the same as in Fig. [ZJ 



• MODSEL: In the block MODSEL, CP violation with completely general CP phases has 
been selected. 

• SMINPUTS_and VCKMIN: The quantities a-£(M z ) m , G F , a s (M z ) m , M z , m b (rnb) m , 
m c (m c ) MS , pole masses of top-quark, electron and muon, and md.u,s{^ GeV) MS are 
given in the SMINPUTS block. In the VCKMIN block, the four parameters for the CKM 
mixing matrix A, A, p, and f] are given. 

• EXTPAR and IMEXTPAR: The blocks for non-minimal parameters have been filled ac- 
cording to Section 2 of the SLHA2 writeup [11] 0. For the input scale M inpu t, we are 
taking Q tb , the scale of the heaviest third-generation squark |20j . which is stored in 
RAUX_H(13)= Ql b . The imaginary parts of the gaugino masses, the trilinear couplings, 
and the \i parameters are listed in the block IMEXTPAR. 

• MASS: In the MASS block, the pole masses of the neutral and charged Higgs bosons 
are given, and the iy-boson mass is calculated from the Fermi constant via M\y = 

^For the Higgs parameters, only the relevant quantities ^(Afi nput ), tan/3(Mi nput ), and the charged Higgs 
pole mass have been included. 
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Figure 8: The LEP and LHC exclusion regions in the MHMAX scenario. 



1/2 



i^g 2 1 '4v / 2G fJ ' . Note that, in the present version, the treel-level masses of the third- 
generation sfermions, neutralinos, and charginos are given. 

HCOUPLINGS, IMHC0UPLINGS: In these blocks, the seven Higgs-self couplings Ai_ 7 cal- 
culated according to Ref. [21] are given from CAUX_H(201— 207) . 



THRESHOLD: In this block, the measures of the threshold corrections to the Yukawa 
couplings h t .b are given from CAUX_H (21 1,212). The quantity A& is also available 
from CAUX_H(10). 

NMIX, IMNMIX; UMIX, IMUMIX, VMIX, IMVMIX; ST0PMIX, IMST0PMIX; SB0TMIX, IMSB0TMIX; 
STAUMIX, IMSTAUMIX: In these blocks, the mixing matrices of neutralinos, charginos, 
sfermions are given. Precisely, 

- (NMIX + / IMNMIX) ja = Nu, with i = 1, 2, 3, 4 and a = (B, W 3 , S}, Sj) , 

- (UMIX + / IMUMIX) to = (C L ) ia with i = 1, 2 and a = (W~, H~) L , 

- (VMIX + / IMVMIX) to = (C R )* with i = 1, 2 and a = (W~, H~) R , 



(ST0PMIX + 7 IMSTDPMIX) ia = 
larly for sbottoms and staus. 



IP 



with i — 1,2 and a = (tL,tR), and simi- 
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Figure 9: The upper frames are for the decay widths of the SM Higgs boson into b quarks 
(upper left) and photons (upper right) as functions of its mass. The lower frames are for 
the branching ratios into W (lower left) and Z (lower right) bosons. In the lower frames, 
the solid (dotted) lines are obtained using 4 (3)-body phase space below the mass thresholds. 
The vertical lines correspond to Mh sm = 125.5 GeV. 



CVHMIX: The block for the neutral Higgs boson mixing has been filled as suggested in 
Ref. HI]: 





( M 


(CVHMIX)^ = 


oL 3 o 




I /. 



\ 



( 1 

10 

— sin /3 cos j3 

V cos/3 sin/3 / 



(20) 



with i = 1, 2, 3 



a 



>i, (/> 2 , a, G°) and a 



h, $2,0,1,0*2) 



We note that a 



sin a\ + cos /3 02 and G° = cos a\ + sin /3 02- 



AU, IMAU, AD, IMAD, AE, IMAE; YU, IMYU, YD, IMYD, YE, IMYE: In these blocks, the 
third-generation (3, 3) components of the trilinear-coupling and the Yukawa-coupling 
matrices are given at the scales Qtb and mf ole , respectively. 

DECAY: The total decay widths and the non- vanishing branching ratios of the neutral 
and charged Higgs bosons are stored in the block from the arrays GAMBRN(101 , 1 , IH) 
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CPX scenario with tan(3=5 and <I> A = <& 3 = 90° 





Figure 10: The four Mercury EDMs and the Radium EDM as functions of the parameters 
p that parameterizes the hierarchy between the first two and third generations as m^ i2 = 
pmx 3 with X = Q,U, D, L, E . The CPX scenario has been taken with tan/3 = 5, $A t6 = 
$ 3 = 90°, and M SU sy = 0.5 TeV: = = = m Ls = rng 3 = M SU sy/ M = 

4 Msusy , \A t , b ,r\ = 2 Msusy , \M 3 \ = 1 TeV; M 2 = 2M X = 100 GeV , M H ± = 300 GeV. 

and GAMBRN(IM,3,IH) with IH= 1 -3 and GAMBRC(51 , 1) and GAMBRCQM, 3), respec- 
tively. The decay width and branching ratios of the top quark are also stored from 
the arrays RAUX_H (50-53) . 

• FOBS: In this block, we show the branching ratios B{b —t X 8 j), B(B S — > fi/j,), 
B(Bd — > tt), the ratio Rbtu, and the CP asymmetry Acp(B — > X s j). These quanti- 
ties are also available from RAUX_H(130, 131) and RAUX_H(134, 135, 136) in different 
normalizations. 

• F0BSBSM: In this block, we show the SUSY contributions to AM Bd ,B s from RAUX_H( 132, 133) . 
For the total AM Bd)Ba , one may add the quantities (B®\ H^ B=2 \B^) S u SY and 

(B° s \ H^ B=2 \B° S ) SVSY given in CAUX_H(150 , 151) to the SM contributions. 

• FDIP0LE: Here, we list the EDMs of Thallium, neutron, Mercury, Deuteron, Radium, 
and muon, as well as the anomalous magnetic moment of muon, (g^ — 2). In the cases 
of the neutron and Mercury EDMs, we present 3 and 4 evaluations, respectively, for 
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the estimation of the theoretical uncertainties. 



• HiggsBoundsInputHiggsCouplingsBosons, HiggsBoundsInputHiggsCouplingsFermi 

In these blocks, the Higgs couplings normalized to the corresponding SM couplings 
are listed for the interface to the program HiggsBounds [15J. 

6 Summary 

Encouraged by the recent observation of a new particle resembling a SM-like Higgs boson 
at the CERN Large Hadron Collider [J], in CPsuperH2.3 we have performed a number 
of updates to the CPsuperH code. In detail, we have improved the computation of CP- 
violating effects on Higgs-boson masses and mixing, by including stau contributions [5] and 
finite radiative effects on the tau-lepton Yukawa coupling. These effects play an important 
role in the decays #1,2,3 — > 77; for large values of tan/3 [TU|22]. We have also implemented 
the LEP limits on the processes e + e~ — > HiZ, HiHj, including the bounds on Hi —> fr 
obtained by CMS with 4.6 fb _1 of LHC data at a centre-of-mass energy of 7 TeV. Finally, 
we have included the decay modes of the neutral Higgs bosons #1,2,3 into a Z boson and 
a photon. These decay modes are expected to become observable as more luminosity is 
accumulated at the LHC. 

We have also incorporated in CPsuperH2 . 3 calculations of the EDMs of Mercury and 
225 Ra, including estimates of the contributions due to Schiff moments. We have presented 
a number of numerical examples and figures for each of our updates, in order to illustrate 
some typical results obtained by CPsuperH2 . 3. To enhance the synergy and compatibility 
of CPsuperH2 . 3 with other codes, we have created an SLHA2 interface in accordance with 
the SUSY Les Houches Accords. For comparison with other works, we include evaluations 
of the principal decay rates and branching ratios of a SM Higgs boson. In the Appendix, 
we exhibit a number of output Tables, in order to highlight the updates with respect to 
the older version of CPsuperH. 

If indeed the new particle discovered by ATLAS and CMS turns out to be a Higgs 
boson, the joint task of theory and experiment will be to establish whether it is compatible 
with the SM or exhibits deviations characteristic of some extension of the SM such as 
the MSSM. The possibility of CP violation should be considered in studying the latter 
possibility, and CPsuperH2 . 3 is a suitable updated tool for this task. 
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Appendix 



Most importantly, in addition to the SLHA2 interface, the array for the SM parameters 
SMPARA_H has been extended to include the SM Higgs mass,M# SM , see Table [TJ Because of 
these improvements, from CPsuperH2 . 2 to CPsuperH2 . 3, the following changes to the files 
run and cpsuperh2 . f are needed: 

• In the run file: the entries SMPARA_H(20) and IFLAG_H(30) added 

• In the cpsuperh2 . f file: 

- REAL* 8 SMPARA_H(19) , SSPARA_H(38) =>• REAL*8 SMPARA_H(20) , SSPARA_H(38) 



DATA NSMIN/19/ 



DATA NSMIN/20/ 



Furthermore, in the updated version, the contents of the array GAMBRN(IM, IWB=2 , IH) for 
the neutral Higgs-boson decays now include the corresponding SM branching ratios. To 
reproduce Figs. |6] and [TJ for example, one simply needs to use the following ratios of 
parameter arrays: 



B(Hi — > 77) 



MSSM 



GAMBRN(IM = 17, IWB = 3, IH = 1) 



B(H l ->■ 77) sm 
B(Hi — > ^7) mssm 



GAMBRN(IM = 
GAMBRNf IM 



17, IWB = 2, IH = 
19, IWB = 3, IH 



B{H X -> Zj) 



SM 



GAMBRN(IM = 19, IWB = 2, IH = 1 



For the calculation of BiHi — y Zj), we closely follow Ref. [27] . 

The amplitude for the decay process Hi — y Z(ki, ei) 7(^2, £2) can be written as 

a 



M 



2irv 



{S?\M Hi ) [h-hel-el-h-elh-el] - Pf y (M Hi ) (e^A*)} 



(A.l) 

where k^ are the momenta of the Z boson and the photon (we note that 2k\ ■ k 2 = 
Mfj. — Mf ), ei j2 are their polarization vectors, and (e^/ci/^) = e^a/^ie^i ^2 • The decay 
width is given by 



T(H t -> Zi) 



aGlM^ w _ M ^ ( Ml 



64tt 4 

1287T 3 U 2 



F?\M Hi )~+ F?>(M Hi 



M 

Ml. 



2 \ 3 



Fi Q \M H y+ F?>(M Hi 



(A.2) 
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The scalar and pseudoscalar form factors are given by 



£ Af+A^l+A w + A H± + £ Aj, 

f=t,b,T f=t,b,T 

I (5) , A ( 5 ) 



E + ^ 

f=t,b,r 



= ■ 



with 



4(0),(5) 



^3 ^ s wQ'f 



A 



(0),(5) 

Ah± 
A 



The loop functions are El 



9D A^rr, 2 J 3 "lyV/ 5,P p (0),(5) 

- M§ cot 9 W F~±^ , 
M| cot 9 w g HiVV F w , 

,2 

— gHiH+H-F H ± , 

M 2 N f c Q f Fj. 



v 



F (5) 

Fh ± 

F f 



C (m 2 )+AC 2 (m 2 ), 
Co(mJ) , 

m~± / x 
2^ £^/(m~ f ,m~±,m~±J 



2 ^ Y.Tzr.9\ m xt^ m xi^ m xi 



M w 



z x +x- k g H, x + x - ' 
v Zxtx~jH l xtx- ' 



Hi { 1 Q „2 



l-2c^) + 2(l-6c 2 w ) 



C 2 (M^)+4(l-4c 2 w/ )C (M^) 



4C 2 (M^ ± ) 
4v 2 



(A.3) 



(A.4) 



(A.5) 



We follow the conventions and notations of CPsuperH [2] for the Higgs couplings to the 
SM and supersymmetric particles, and the relevant Z-boson interactions are given by the 
following Lagrangian terms: 



Z-sfermion-sfermion 



C 



z.ff 



where gz = e/(swcw) and 



with I^' v = +1/2 and I : 



d,e 



HutM - Q f s 2 w 8 l3 



(A.6) 



(A.7) 



4/2. 



**For the functions of Co,2(w 2 ), f (mi, 1712,1712), g(rrii, 1712,1712), and 02(1711,1712,1712), we refer to [27] . 
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Z-chargino-chargino [28] 



C 



Zx+x~ 



9zX l T \v Zxt 



(A.8) 



where 



J zxtx] 



a zxfx~ 



4 

1 r 



(Cl) i2 (C l )* 2 + (C R ) a (C R ) 

(c L ) l2 (c L )% - (c R y i2 (c R ) 



.ft 



J'2 



(A.9) 



For completeness, we recall that the Z-boson couplings to the quarks and leptons are 
given by 



C 



Zff 



9z 



(A.10) 



with v Z f f = 4/2 - Q f s 2 w and a z j f = l(/2. 

B Tables updated 

In this Appendix, we provide output Tables extended to include the updates implemented 
since the appearance of CPsuperH2 . 0. 

Table 1: The contents of the extended SMPARA_H(IP) . 



IP 


Parameter 


IP 


Parameter 


IP 


Parameter 


IP 


Parameter 


1 

2 
3 
4 
5 


a s {M z ) 

M z 
sin 2 9w 

m e 


6 
7 
8 
9 
10 


m T 
m d (m v t ole ) 

m s (mr' ) 
m b (m\ ole ) 


11 
12 
13 
14 
15 


m u {m p t ole ) 
m c (m£ ole ) 

pole 

T w 


16 
17 
18 
19 
20 


A 
A 

P 
fj 

Mh sm 
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Table 2: The contents of the extended SSPARA_H(IP) 



IP 


Parameter 


IP 


Parameter 


ip 


Parameter 


IP 


Parameter 


1 


tan/3 


11 


m Q, 


21 




31 


14.1 


2 


MP°i c 


12 


m u 3 


22 


PQ 


32 


$A U 


3 


H 


13 




23 


Pu 


33 


141 


4 


% 

m\ 


14 


m u 


24 


Pd 


34 


$A C 


5 


15 


m E 3 


25 


PL 


35 


141 


6 




16 


|4| 


26 


Pe 


36 




7 


|M 2 | 


17 


®A t 


27 


141 


37 


141 


8 




18 


|4| 


28 




38 


®A S 


9 


|M 3 | 


19 




29 


14*1 


39 




10 


$ 3 


20 


|4| 


30 


®A» 


40 
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Table 3: The contents of the array RAUX_H. In RAUX_H(22) and RAUX_H(23), the notation 
h® is for the Yukawa couplings without including the threshold corrections. The notations 
are explained in Refs. U\ \B, IWj \2( \ [2^2hJ . For EDMs and the magnetic dipole moment 
(MDM) of the muon, some following slots of the array are allocated for the constituent 
contributions, see Refs. IWj , 125^ and I2b^ for details. 



pole 

ml 

/ pole\ 



RAUXJKl) 
RAUX_H(2) 
RAUX_H(3) 
RAUX_H(4) 
RAUX_H(5) 
RAUX_H(6) 



RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 
RAUX. 



(10 
(11 
(12 
(13 
(14 
(15 
(16 
(17 
(18 
(19 
(20 
(21 
(22 
(23 
(24 
(25 
(26 
(27 
(28 
(29 
(30 
(31 
(32 
(33 
(34 
(35 
(36 



a 



:« 0le ) 



pole 



m 
m 

a Am 



mP° le ) 

pole") 



(mP° le ) 



Mff or Mf ± 



Qi 
Ql 



/ polc\ 

v\[m% ) 
Vi(Q t ) 
vi(Qb) 
vi(Qtb) 

v 2 {Qt) 
v 2 (Qb) 

V2(Qtb) 

/^« ole )l 
/^K ole )l 

h t (mr lc )\ 
ht(Qt)\ 
h t (Qtb)\ 

M m r lc )i 

h(Q b )\ 
h{Qtb)\ 
M\ 

$lefl H+H ~ (M 
\jv\mr lc )/2 
A 4 (m polc ) 
Ai (mr lc ) 

A 2 K ole ) 

A34W 0le ) 



pole 2 ^ 



RAUX_H(101) 



RAUX 
RAUX 
RAUX 
RAUX 
RAUX 
RAUX 
RAUX 



.H(130) 
.H(131) 
.H(132) 
.H(133) 
.H(134) 
.H(135) 
.H(136) 



B{B s ^fifi)xl0 7 
B(B d ->TT) xlO 7 

Rbtu 

B(B^X sl )xW 



RAUX_H(200) (df/e) cm 

RAUX_H(210) (rff/e) cm 

RAUX_H(220) (df/e) cm 

RAUX_H(230) (df/e) cm 

RAUX_H(240) d% cm 

RAUX_H(250) d c d cm 

RAUX_H(260) d G cm/GeV 

RAUX_H(270) C s cm/GeV 

RAUX_H(271) C P cm/GeV 

RAUX_H(272) C' P cm/GeV 



RAUX 
RAUX 
RAUX 
RAUX 
RAUX 
RAUX 
RAUX 
RAUX 
RAUX 
RAUX 
RAUX 



JK280) 
JK281) 
.H(282) 
.H(283) 
.H(284) 
.H(285) 
.H(286) 
.H(287) 
.H(288) 
.H(289) 
.H(290) 



C de /m d cm/GeV 2 
C se /m s cm/GeV 2 
C e d/m d cm/GeV 2 
C es /m s cm/GeV 2 
C eb /m b cm/GeV 2 
C ec /m c cm/GeV 2 
C et /m t cm/GeV 2 
C dd /m d cm/GeV 2 
C sd /m s cm/GeV 2 
C bd /m b cm/GeV 2 
C db /m b cm/GeV 2 



RAUX_H(300) d T1 ecm 

RAUX_H(310) d n CQM ecm 

RAUX_H(320) d^ M ecm 

RAUX_H(330) (f n QCD ecm 

RAUX_H(340) d Hg ecm 

RAUX_H(350) d D ecm 

RAUX_H(360) {dje) cm 

RAUX_H(380) (a M )susY 



RAUX_H(50) 
RAUX_H(51) 
RAUX_H(52) 
RAUX_H(53) 



T(t W+b) 
T(t -> H+b) 
B{t -> W + b) 
B(t ->■ H+b) 



RAUX_H(400) d c e cm 



RAUX_H(410) 
RAUX_H(411) 
RAUX_H(412) 
RAUX_H(413) 
RAUX_H(414) 



o?Hg ecm 
(ijj g ecm 
djj g ecm 



e cm 



7111 



dft g ecm 



RAUX_H(420) d Ra ecm 



RAUX_H(430) ILEP 



RAUX_H(440) ILHC7f 6 ^ rr 



QSr 



Table 4: The contents of the array RAUX_H, continued. 



RAUX_H 


[501) 


m b (M Hl ) 


RAUX_H(600) 


Mh sm 










RAUX_H 


;502) 


m t (M Hl ) 


RAUX_H(601) 




-> ee) 


RAUX_H(651) 


B(H SM 


-> ee) 


RAUX_H 


;503) 


m c {M Hl ) 


RAUX_H(602) 


r(-^sM 


— > /i/i) 


RAUX_H(652) 


B(H S m 


—> /i/i) 


RAUX_H 


:504) 


m b (M H2 ) 


RAUX_H(603) 




— > tt) 


RAUX_H(653) 


B(H SM 


-»■ tt) 


RAUX_H 


;505) 


m t (M H2 ) 


RAUX_H(604) 




— > old) 


RAUX_H(654) 


B(H SM 


— > dd) 


RAUX_H 


;506) 


m c (M H2 ) 


RAUX_H(605) 


r(-^sM 


-> ss) 


RAUX_H(655) 


B(H S m 


-> ss) 


RAUX_H 


:507) 


m b (M Hz ) 


RAUX_H(606) 


f(h sm 


^bh) 


RAUX_H(656) 


B(H SM 


^bb) 


RAUX_H 


;508) 


m t (M H3 ) 


RAUX_H(607) 


r(-^sM 


UU) 


RAUX_H(657) 


B(H S m 


—> uu) 


RAUX_H 


;509) 


m c (M H3 ) 


RAUX_H(608) 


r(-ff~sM 


-> cc) 


RAUX_H(658) 


B(H SM 


-> cc) 








RAUX_H(609) 


r(-ffsM 


-> tt) 


RAUX_H(659) 


B(H SM 


-»■ tt) 


RAUX_H 


:5ii) 


M Hl (f) 


RAUX_H(610) 


r(-ffsM 


->■ WW) 


RAUX_H(660) 


B(H$m 


->■ WW) 


RAUX_H 


:5i2) 


M H2 (f) 


RAUX_H(611) 


r(-ffsM 


->■ ZZ) 


RAUX_H(661) 


B(H S m 


->■ ZZ) 


RAUX_H 


:5i3) 


M H3 (f) 




















RAUX_H(617) 


r(i?sM 


-> 77) 


RAUX_H(667) 


B(H SM 


->■ 77) 


RAUX_H 


;520) 


OUT) 


RAUX_H(618) 


r(£?sM 




RAUX_H(668) 


B(H SM 




RAUX_H 


[521) 


OUT) 














RAUX_H 


[522) 


OUT) 


RAUX_H(650) 


pSM 
1 tot 




RAUX_H(700) 


e>SM _ 
^tot — 


1 


RAUX_H 


[523) 


OUT) 














RAUX_H 


[524) 


OUT) 














RAUX_H 


;525) 


OUT) 














RAUX_H 


[526) 


Oal{f) 














RAUX_H 


[527) 


o*(f) 














RAUX_H 


;528) 


Oas(f) 
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Table 5: The contents of the array CAUX_H. The notations are explained in Refs. IHHF 
andfEfl. For CAUX_H( 10, 12), h q = h° q /(l + A g tan/3). 



CAUX_H(1) 


ht/\ht\ 


CAUX. 


_H( 


:i30) 






CAUX_H(2) 


h/\h\ 


CAUX. 


_H( 


:i3D 


p?(V§) 




CAUX_H(10) 


A fe 


CAUX. 


_H< 


:i32) 






CAUX_H(11) 


A s 


CAUX. 


_H< 


:i33) 






CAUX_H(12) 


h s (m\ olc ) 


CAUX. 


_H( 


;i34) 


Sj{yfl) 




CAUX_H(13) 


h T 


CAUX. 


_H( 


:i35) 


P 3 7 (v / I) 




CAUX_H(100) 














CAUX_H(101) 




CAUX. 


.H( 


:i4o) 


Sf(VB) 




CAUX_H(102) 




CAUX. 


.H( 


[141) 


pf(Vs) 




CAUX_H(103) 




CAUX. 


.H( 


[142) 






CAUX_H(104) 


D 2,l (*) 


CAUX. 


_H( 


[143) 


P£(V§) 




IS TT"\7" TT /" A /"\ 1 — \ 

CAUX_H(105) 




CAUX. 


_H( 


[144) 


Si(V2) 




H ATTY TJ ( -1 C\G.\ 


r>H°f$\ 

JJ 2,3 \ S ) 


PATTY 


.n \ 


. 14o J 


Pi(V§) 




CAUX_H(107) 


2,4 V / 












CAUX_H(108) 




CAUX. 


Si 


[150) 


m\H*k B = 


2 -B°)susy 


CAUX_H(109) 




CAUX. 


.H( 


[151) 




2 1 ro\ 


CAUX_H(110) 














CAUXJH(lll) 




CAUX. 


.H( 


[160) 


c (0)SM (m pol C) 


CAUX_H(112) 


4,1 V / 


CAUX. 


.H( 


[161) 


Cf SM (mr lc ) 


CAUX_H(113) 


£>g(S) 


CAUX. 


_H( 


[162) 


C (0)SM (m pol e) 


CAUX_H(114) 




CAUX. 


.H( 


[163) 


n (Q) SM+^± 


«° lC ) 


CAUX_H(115) 




CAUX. 


.H( 


[164) 


^,(0) SM+J/± 

o 7 


«° 1C ) 


CAUX_H(116) 




CAUX. 


_H< 


[165) 


n {0)SM+H ± 


«° 1C ) 


CAUX_H(117) 




CAUX. 


Si 


[166) 


c< (0)SM + ^ + x± (m pole ) 


CAUX_H(118) 




CAUX. 


_H( 


[167) 


C? )SM+H±+ji± (mF le ) 


CAUX_H(119) 




CAUX. 


.H( 


[168) 




+ ^(mr le ) 






CAUX. 


Si 


[169) 


Cs(mi° le ) 


1/GeV 






CAUX. 


Si 


[170) 




1/GeV 






CAUX. 


Si 


[171) 


C 10 (ml° le ) 
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Table 6: The contents of the array CAUX_H, continued. 



CAUX_H(201) 


Ai 






CAUX_H(202) 


A 2 






CAUX_H(203) 


A 3 






CAUX_H(204) 


A 4 






CAUX_H(205) 


A 5 






CAUX_H(206) 


A 6 






CAUX_H(207) 


Av 






CAUX_H(211) 


A 6 tan/3/(l - 


f- A fe tan/3) 




CAUX_H(212) 


A t cot/3/(l H 


- A t cot 0) 




CAUX_H(221) 


S 9 SM (M Hl ) 






CAUX_H(222) 


Slu(M H2 ) 






CAUXJK223) 


S 9 SU (M H3 ) 






CAUX_H(231) 


SUM Hl ) 






CAUX_H(232) 








CAUX_H(233) 


SUMh 3 ) 
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Table 7: The updated contents of the array GAMBRNQM, IWB=1 , IH) for the decays of the 
neutral Higgs bosons H m . The entry GAMBRNQM, IWB=1 , IH) is for the decay width of the 
decay mode IM in GeV. And the entries GAMBRNQM, IWB=2 , IH) and GAMBRNQM, IWB=3, IH) 

are the corresponding SM and the full SUSY branching ratios, respectively. 



TM 
±11 




TM 
XI 1 




TM 
J. 11 


1 iDfQV \/l/~\rid 


1 

± 


11 IH ^ t t 


1 1 

X J. 


J.Z jfj r Zj Zj 






2 


Htu — V UU 
J.H ' ft* 


12 


H ttj — ^ H i /j 

In ' 1 






3 


H m -> rf 


13 


-f^iH H 2 Z 






4 


H m — > dd 


14 


Hm — > H 1 H 1 






5 


H m — > ss 


15 


Hm — > H 1 H 2 






6 


H m bb 


16 


Hm — > H 2 H 2 






7 


H m — > uu 


17 


H m -> 77 






8 


H m — > cc 


18 


Hm^ 99 






9 


H m ti 


19 


H 1E ^Z 1 






10 


H m -> WW 






ISMN 


GAMSM 


IM 


Decay Mode 


IM 


Decay Mode 


IM 


Decay Mode 


ISMN+1 


^ih xix°i 


ISMN+11 


#ih -> xtxT 


ISMN+21 


H IE -> &2&1 


ISMN+2 


H m xixi 


ISMN+1 2 


Hm -)> Xi"X2 


ISMN+22 


Fih ->■ b 2 b 2 


ISMN+3 


Hm -> 


ISMN+1 3 


#IH X2~XT 


ISMN+23 


H m -»■ 7?ri 


ISMN+4 


#ih Xi£! 


ISMN+1 4 


#IH x4X2 


ISMN+24 


#IH -> 7?f 2 


ISMN+5 


#IH xlxl 


ISMN+1 5 




ISMN+25 


#ih -> r|ri 


ISMN+6 


#ih xlxl 


ISMN+1 6 


-fflH — > t\t 2 


ISMN+26 


#IH -> 72^2 


ISMN+7 


Hm x1x2 


ISMN+1 7 


Hm t 2 ti 


ISMN+27 


#IH -> K V t 


ISMN+8 


#ih -> xlx! 


ISMN+1 8 


Hm t^t 2 






ISMN+9 


#ih X3X4 


ISMN+1 9 




ISMN+ISUSYN 


GAMSUSY 


ISMN+10 


#IH X°X° 
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